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ABSTRACT 

We present high-resolution, near-infrared imaging of the circumnuclear region of 
the barred spiral galaxy M100 (=NGC 4321), accompanied by near-infrared spec- 
troscopy. We identify a total of 43 distinct regions in the iC-band image, and determine 
magnitudes and colours for 41 of them. By comparison with other near-infrared maps 
we also derive colour excesses and X-band extinctions for the knots. Combining the 
imaging and spectroscopic results, we conclude that the knots are the result of bursts 
of star formation within the last 15-25 Myr. We discuss the implications of these new 
results for our dynamical and evolutionary understanding of this galaxy. 

Key words: galaxies: evolution - galaxies: individual (M100, NGC 4321) - galaxies: 
kinematics and dynamics - galaxies: structure - infrared: galaxies. 



1 INTRODUCTION 

, The circumnuclear regions (CNRs) of barred spiral galax- 

■ ies, with or without an active galactic nucleus, often show 
significant star forming activity, frequently organized into 
ring-like structures (e.g., Morgan 1958; Sersic & Pastoriza 
1967; Pogge 1989a, b; Buta & Combes 1996). The massive 
star formation (SF) is usually observed by imaging in the 

■ Ha emission line, and is understood to occur as a result 
of orbit crowding and gas accumulation between a pair of 
dynamical Inner Lindblad Resonances (ILRs) in the disc. Al- 
though the SF in the circumnuclear 'rings' can be explained 
by gravitational instabilities there (Elmegreen 1994), spiral 
structure is almost always present, either in the SF itself 
(e.g., M100: Pogge 1989a; Cepa & Beckman 1990; Knapen 
et al. 1995a, hereafter Paper I), or in dust lanes associated 
with the CNR (e.g., NGC 5248: Laine et al. 1998), indicat- 
ing that the SF is organized, e.g., by density waves. Near- 
infrared (NIR) imaging of CNRs gives important clues on 
the distribution of obscuring dust (Paper I), and can also 
be used to derive the distribution of stellar mass, which in 
turn is used as input for dynamical modelling (e.g., Quillen, 
Frogel & Gonzalez 1994; Wada, Sakamoto & Minezaki 1998; 
Garci'a-Burillo et al. 1998). The NIR light distribution is 
not used as input to models by e.g., Knapen et al. (1995b; 
hereafter Paper II), Wada et al. (1998), and Yuan & Kuo 
(1998), but these authors do compare the output of their 
models with the observed NIR morphology. In Paper II, we 
estimated the value of the Tf-band mass-to-luminosity ratio 
(M/L K ) in the CNR of M100, concluding that M/L K is not 



constant, and is strongly affected by SF. Wada et al. (1998) 
subsequently confirmed this from a comparison of their two 
models. 

The main interpretational issue in this respect is what 
fraction of the NIR emission traces the stellar mass distri- 
bution, and what fraction is from dynamically young stars? 
In Paper II this issue was discussed qualitatively in detail, 
concluding that in the CNR of M100 the latter fraction was 
significant. Whereas the K light is distributed smoothly as 
compared to blue light or Ha in M100 (cf. e.g. NGC 4303: 
Elmegreen et al. 1997), Knapen et al. (1995b) drew atten- 
tion to two A"-band 'hot spots', named Kl and K2, which 
were identified as starburst regions at the loci where a pair 
of trailing and leading shocks interact in the vicinity of the 
inner ILR. Wozniak et al. (1998) subsequently showed that 
at 6.75 and 15 pm, Kl and K2 emit similarly, as expected 
in such a scenario. Knapen et al. also reported the detection 
of 'leading arms' in K, in the inner bar region, predicted 
theoretically. 

In this paper, we quantitatively confirm previous sus- 
picions concerning the stellar populations in the CNR of 
M100, based on an analysis of high-resolution NIR imag- 
ing of the CNR, and of NIR spectroscopy of Kl, K2, and 
the nucleus. After describing the observations (Sect. ^ and 
results (Sect. |^), we discuss the ages of the compact knots 
discovered from the NIR imaging, and the stellar popula- 
tions within them. In Sect. we also discuss how the new 
results ht in with previous work on the CNRs in M100 and 
other galaxies. We conclude briefly in Sect. 0. 



© 1998 RAS 



2 S. Ryder and J. Knapen 



2 OBSERVATIONS 

2.1 High-resolution imaging 

The nuclear region of M100 was imaged with IRCAM3 on 
the upgraded 3.8-m United Kingdom Infrared Telescope 
( UKIRT) and a (warm) magnifier, which yielded a pixel 
scale of 0.140 arcsec and a 36 arcsec square field of view. 
Images in J, H, and K were obtained on the night of 
1998 February 14 UT; bracketing observations of PSF stars 
yielded seeing FWHM of 0.34 arcsec at K and 0.39 arcsec 
at J. For each filter, a 3 x 3 mosaic of overlapping frames 
spaced 4 arcsec apart was obtained on M100, together with 
a matching mosaic on the sky 10 arcmin away. The integra- 
tion time per frame was 60 sec. After alignment, the frames 
in each mosaic were combined, excluding known bad pixels, 
and the sky mosaics were used to fiatfield the object mosaics. 
Owing to their small field coverage, the image mosaics had 
to be flux calibrated using aperture photometry from the 
literature (Devereux 1989; Recillas-Cruz et al. 1991). 



2.2 Near-infrared spectroscopy 

Longslit H and K band spectra were obtained separately on 
the night of 1998 February 12 UT using CGS4 on UKIRT. A 
slit width of 1.2 arcsec was used with the 40 line mm -1 grat- 
ing in first order, yielding resolving powers (after smooth- 
ing) at H and K of 285 and 360, respectively. The slit posi- 
tion angle was set to 114°, allowing us to capture both the 
A'-band hotspots Kl and K2 identified in Paper I, as well 
as the nucleus in the same observation. Since the infrared- 
emitting central region of M100 fills a sizable fraction of 
the 90 arcsec slit length, sky spectra were taken with the 
object positioned well off the slit, rather than simply dis- 
placed along the slit. Total on-source integration times were 
1920 sec and 2400 sec at H and K, respectively. Observa- 
tions of the A4 V star BS 4632 and the F0 V star BS 4531 
were used to divide out atmospheric absorption features in 
the M100 spectra, after first having their own Brackett and 
Paschen lines removed. Preliminary image processing (sub- 
traction of object-sky pairs and flatfielding) was performed 
using the on-line reduction software CGS4dr, while object 
extraction, wavelength calibration, and division by the stan- 
dard were done with irafQ. 

3 RESULTS 
3.1 Photometry 

Figure ^ shows the reduced if-band image. A total of 
43 compact 'knots' were identified by a comparison of this 
image with an unsharp-masked version (cf. Mighell & Rich 
1995). The improvement in delivered image quality as a re- 
sult of the UKIRT Upgrades Program (which includes active 
control of the primary mirror figure and telescope focus, a 
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tip-tilt Fast Guider, and dome ventilation louvres; Hawar- 
den et al. 1998) can be clearly seen by comparing this im- 
age with the pre-upgrade UKIRT image of Paper I (their 
fig. 2). What were previously just broad maxima in the K- 
band light have now been resolved into a host of individual 
regions. For example, both the maxima Kl and K2 identi- 
fied in Paper I at 0.8 arcsec resolution turn out to consist of 
3 or 4 sub-clumps at 0.4 arcsec resolution. The detection of 
leading arms within the SF 'ring' (Paper I) is confirmed on 
the new images. 

Photometry of each of these knots in all 3 bands was 
performed using the CCDCAP digital circular aperture pho- 
tometry code developed by Mighell to analyze HST WF/PC 
and WFPC2 observations (Mighell 1997, and references 
therein). The local (sky + unresolved bulge light) back- 
ground was determined from a robust estimate of the mean 
intensity value of all pixels between 15 and 20 pixels of the 
aperture centre. To minimise the effects of crowding, magni- 
tudes were determined by summing the flux within a radius 
of 5 pixels (0.7 arcsec), then corrected to an aperture radius 
of 15 pixels based on similar measurements of the PSF stars. 
The positions, Tf-magnitudes, and NIR colours (uncorrected 
for reddening), colour excesses and inferred Tf-band extinc- 
tions (see below) for the 41 knots that could be reliably 
measured are tabulated in Table |l[ 

As is readily apparent from a comparison of Fig. |l| with 
fig. 4 of Paper I, many of these knots could be severely red- 
dened by their proximity to the dust lanes. We have made 
a first-order correction for this in the following way. Us- 
ing the isochrone synthesis models of Bruzual & Chariot 
(1993), we find that from an initial (J — K) colour of 1.48 at 
20 Myr, a given population will become bluer by ~ 0.1 mag 
at 1 Gyr, then finish up redder by about the same amount 
after 10 Gyr. We have therefore used the (/ — K) map from 
Paper I in order to calculate the colour excess (relative to 
1.48) of the unresolved stellar background around each of 
the knots, together with the interstellar extinction curve of 
Rieke & Lebofsky (1985), to derive appropriate corrections 
for reddening and extinction at the location of each of the 
knots, as given in Table |l| 

The distribution of knot NIR colours, after correcting 
for reddening as above, is shown in Fig. |2| The bulk of the 
knots cluster around a point slightly bluer than the mean 
colours for an old stellar population, i.e., 0.71-0.76 in (J — 
H), and 0.25-0.30 in (H - K) (Griersmith, Hyland & Jones 
1982). Some of the knots are redder in (J — H) and bluer 
in (H — K) than can be accounted for by their photometric 
errors, or by any amount of internal reddening, warm dust, 
or scattered light (see, e.g., Israel et al. 1998). These are 
all diffuse, outlying regions. The three bluest knots (33, 34, 
and 42 in Table hj) are all located along the inner edge of 
the southwest dust lane. 

The A"-band Luminosity Function (LF), corrected for 
extinction using the Ak values in Table jj], is presented in 
Fig. |. A distance to M100 of 16.1 Mpc (Ferrarese et al. 1996) 
has been assumed. The LF peaks at Mk ~ —14, and has a 
steep bright-end cutoff compared with the faint-end tail. We 
have overlaid on Fig. ^ the characteristic LF for extragalactic 
globular cluster systems, a Gaussian having (Mv) = —7.1 
and a(Mv) ~ 1.3 (Harris 1991). It is possible to fit the 
peak and the faint-end, provided (V — K) — 6.9 (which 
would require the stellar content to be totally dominated by 
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dwarfs or giants of type M6 or later), but there is then a 
clear deficiency of bright objects if, indeed, these knots are 
globular clusters. 

3.2 Spectroscopy 

Continuum fitting and equivalent width measurements were 
carried out using the SPLOT task in iraf. The results are 
given in Table ^. The l-cr measurement errors are ~ 0.1 nm 
for the emission lines. In the case where no line feature was 
apparent at the expected (redshifted by 1571 km s _1 ) wave- 
length, the nearest positive peak was measured and used as 
an upper limit. For the CO measurements, we have used 
the wavelength range W2 of Puxley, Doyon & Ward (1997; 
hereafter PDW), equivalent to that used by Kleinmann & 
Hall (1986) and by Origlia, Moorwood & Oliva (1993). In 
the last column of Table |2| we give the spectroscopic CO 
index (CO ap ) as defined by Doyon, Joseph & Wright (1994; 
hereafter DJW). The use of a broader wavelength window 
here is important, since our resolution is at least a factor of 
2 poorer than results published elsewhere. Kl and K2 have 
quite similar CO line strengths, and both are stronger than 
the nucleus. 

The values observed (~ 1 nm) are similar to those in 
M83 (PDW), and galaxies in general (Oliva et al. 1995). We 
note in passing that according to the diagnostic diagram of 
Moorwood & Oliva (1988), the ratios of H 2 1-0 S(l) to Br7, 
together with the upper limits on [Fen], would tend to sug- 
gest that the nucleus of M100 has some Seyfert 2 activity, 
while Kl and K2 would be classed as 'Composite'-like (i.e., 
their optical spectra would have indications of Seyfert and 
Hn activity). However, we feel that many of the line detec- 
tions are too marginal for the ratios to have much meaning. 



4 DISCUSSION 

4.1 Old stellar clusters or recent massive star 
formation? 

Having resolved these compact, NIR emission knots, we now 
attempt to draw some inferences about their nature. In par- 
ticular, are these knots globular clusters, dominated by an 
old stellar population? Or are they much more recent, per- 
haps formed in a single burst, or triggered sequentially by 
some dynamical process? 

On the basis of the photometry alone, we might be 
tempted to conclude that these knots are a few Gyr old. 
They are only slightly bluer than the characteristic colours 
of an old stellar population. The A"-band LF can be approx- 
imated by the Gaussian globular cluster LF, aside from a 
deficiency of bright objects. It is the addition of the spec- 
troscopy, however, that provides the decisive argument in 
favour of these knots being the result of recent star forma- 
tion. The mere presence of Br7 in emission, coupled with 
the deep CO absorption features attributable to young su- 
pergiants, argues strongly that (at least for the knots asso- 
ciated with Kl and K2), massive star formation has taken 
place within the last 100 Myr or so. 

We can use our knowledge of the colours, together with 
the starburst models of Leitherer & Heckman (1995; here- 
after LH95) to place some constraints on the ages of these 



knots. To begin with, we have tracked the time evolution of 
(J — H) and (H — K) in the LH95 models for an instanta- 
neous burst, and for continuous star formation, using solar 
metallicity, a Salpeter-type IMF, and a mass range of 1— 
100 M©. The end-points of these two models after 300 Myr 
of evolution are shown as the large diamond and cross, re- 
spectively, in Fig. |^. It can clearly be seen that the knot 
colours on the whole are more consistent with a recent burst 
of star formation, than with either continuous or ancient 
star formation. However, as Satyapal et al. (1997) point out, 
extinction-corrected colours alone do not provide good con- 
straints on stellar population models. 

A more precise determination of the age of the burst 
requires spectroscopic information, which is only available 
for the regions Kl and K2 (and in fact, the narrow slit en- 
compasses knots 29 and 19 only). These two regions have 
quite similar colours, emission line equivalent widths, and 
CO spectroscopic indices. Their similarity at mid-IR wave- 
lengths has been noted by Wozniak et al. (1998), and in- 
deed, the main difference between them seems to be the 
degree of obscuration by dust (Paper II). Using Br7 equiv- 
alent widths only, the instantaneous burst models of LH95 
imply an age of 8-13 Myr for these two regions. Even with 
a steeper (Miller-Scalo type) IMF, the LH95 models with 
continuous star formation require > 300 Myr for the Br7 
emission to drop this low. 

Devost & Origlia (1998) have combined the LH95 mod- 
els with equivalent width calculations, isochrones, and evo- 
lutionary tracks from the literature to forecast the time evo- 
lution of the CO (6-3) 1.62 ^m and CO (2-0) 2.29 ^m equiv- 
alent widths. Using their fig. 1, we derive ages of either 5 or 
15 Myr, though the weakness of Br7 in emission makes the 
latter more likely. The DJW models predict the time evolu- 
tion of both Br7 and CO, and by plotting CO sp versus Br7 
line strengths in Fig. ^| (adapted from fig. 3 of PDW), we 
find that both Kl and K2 lie close to the locus traced by 
a starburst evolutionary model having an exponential decay 
time-scale of 5 Myr. The ages of Kl and K2 implied by these 
models are between 15 and 25 Myr. In Paper I, Knapen et al. 
reached similar conclusions regarding the ages of the knots 
Kl and K2 on the basis of (UVK) colours and the models 
of Chariot & Bruzual (1991). 

From modelling and observations, we know that the 
CNR hosts two small-scale but relatively strong density wave 
spiral arms (Paper II; Knapen 1996; Knapen et al., in prepa- 
ration). Starbursts can occur in the spiral armlets, induced 
by the concentration and compression of abundant gaseous 
material by the density wave amplified by the gas self- 
gravity, and show up first as H n regions, and after ~ 5 Myr, 
start emitting strongly in K (LH95; Knapen 1996). The dis- 
tribution of SF regions in the CNR is compatible with and 
fully explained by the interaction of disk gas with the non- 
axisymmetric (barred) potential in M100 (Paper II). Our 
results also confirm that the individual knots of SF in the 
core of M100 have been created in short 15-25 Myr bursts 
by the global density wave pattern. 

4.2 Stellar populations in the knots 

The mean de-reddened colours of the 41 measurable knots 
((./ - H) = 0.74, (H - K) = 0.14) are consistent with the 
stellar population being dominated by giants of type K3-4, 
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or MO supergiants [dwarf stars being too blue in (J — H)]. 
Similarly, the CO sp values for Kl and K2 in Table ti are also 
typical of MO supergiants, or of M4 giants (again, dwarfs do 
not have CO sp this high; DJW). As Oliva et al. (1995) and 
PDW have shown, even the use of CO spectral features at 
1.59, 1.62, and 2.29 /im is still not sufficient to distinguish 
giants from supergiants. The LH95 and DJW models both 
predict the supergiant production rate to peak at ~ 15 Myr 
after the onset of the burst, so a supergiant-dominated pop- 
ulation is perhaps more likely for this reason. In any case, 
we note that a single knot with Mk — — 14 would require up 
to 25000 K4 giants, or as few as just 70 MO supergiants, to 
account for the observed luminosity. Emission lines are not 
expected to contribute much to the knot luminosities, as a 
15 Myr old cluster will have a main sequence turn-off near 
spectral type BO (Whitmore & Schweizer 1995), and such 
stars can no longer ionise the gas around them. 

A difference in the stellar populations of the knots rel- 
ative to the unresolved nuclear light also has implications 
for the dynamical analyses. Our results confirm the varia- 
tions in M/L K over the CNR in M100 derived in Paper II 
and quantified by the dynamical modelling of Wada et al. 
(1998). 

By comparing the K-h&nd luminosities of the LH95 
models 20 Myr after a burst with Fig. ^, we estimate that be- 
tween 0.1 and 4.0 x 10 4 Mq of gas went into forming each of 
these knots, and that the total mass of gas involved in form- 
ing these knots was ~ 6 x 10 6 M© . Since most of these knots 
have been resolved by our observations (FWHM = 4 — 7 pix- 
els), we are able to deconvolve the images, and derive phys- 
ical sizes for the knots of between 30 and 70 pc. 

Finally, we note that these circumnuclear knots in M100 
are somewhat different from the class of 'super star clusters' 
identified in HST images of nearby starburst galaxies (e.g., 
O'Connell, Gallagher & Hunter 1994; Meurer et al. 1995), 
and which have been postulated to be the progenitors of 
globular clusters. The knots are not only larger in size, but 
also much redder than the 'super star clusters', even though 
both types of object are thought to be ~ 15 Myr old. Since 
the nucleus of M100 is not currently undergoing an intense 
starburst, these differences could reflect the different types 
of environment in which they formed. 



5 CONCLUSIONS 

Recent improvements in the delivered image quality at 
UKIRT have allowed us to resolve a number of compact cir- 
cumnuclear 'knots' in M100. NIR colours and luminosities, 
corrected for extinction, are consistent with a recent burst 
(or bursts) of massive star formation as the origin of these 
knots, but the possibility of an older population of objects 
(e.g., globular clusters) cannot be ruled out using this data 
alone. The addition of NIR spectroscopy for two of these 
knots either side of the nucleus provides conclusive evidence 
of recent star formation. By comparing Br7 and CO 2.29 /xm 
line equivalent widths to various starburst models in the lit- 
erature, we infer ages for these knots of between 15 and 
25 Myr, and a stellar population most likely dominated by 
M0 supergiants. Future NIR spectroscopy for the remainder 
of the knots should enable us to generalize this conclusion 
and refine our interpretation of the formation and evolution 



of the SF regions in the presence of the small-scale spiral 
density wave. 
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Table 1. Derived Parameters for Compact Knots in M100 



No. 
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"Offset relative to if-band nuclear position. 

^Insufficient signal-to-noise for centering and photometry. 



Table 2. Line Equivalent Widths in Selected Regions of M100 



Region 


-EW([Fen]) 
(nm) 


-EW(H 2 1-0 S(l)) 
(nm) 


-EW(Br 7 ) 
(nm) 


-EW(H 2 2-1 S(l)) 
(nm) 


EW(CO 2-0) 
(nm) 


CO sp 


Nucleus 


< 0.1 


0.3 


< 0.3 


< 0.2 


0.88 


0.24 ± 0.01 


Kl 


< 0.07 


0.3 


0.4 


< 0.2 


1.01 


0.32 ± 0.02 


K2 


0.2 


0.2 


0.3 


0.1 


0.93 


0.30 ± 0.02 
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Figure 1. /('-band image of the nuclear region of M100, taken with IRCAM3 and a 2x magnifier on UKIRT . Offsets are given relative 
to the isf-band nucleus position. 
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Figure 2. Plot of knot NIR colours, after first-order corrections for reddening as described in the text. The average photometric error 
bar from Table [lj is shown, as is the reddening vector for a residual visual extinction of 2 mag. The open circle marks the 'typical' colours 
of an old stellar population, the open diamond marks the colours of the Leitherer & Heckman (1995) model 300 Myr after a burst of star 
formation, and the hollow cross marks the colours of the Leitherer & Heckman model after 300 Myr of continuous star formation. Filled 
triangles are used for knots embedded within the dust lanes, while filled squares and circles indicate knots adjacent to, or well removed 
from the dust lanes, respectively. 



© 1998 RAS, MNRAS 000, §-10 



Circumnuclear star formation in Ml 00 9 





Figure 3. i^-band Luminosity Function for the 41 measurable objects in Fig. gj, after corrections for extinction. We have superimposed 
on this LF a Gaussian with (Ma) = —13.95, representative of globular cluster systems in nearby galaxies, but there seems to be a 
deficiency of at least 6 bright clusters in this case. 
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Figure 4. Observed CO spectroscopic indices and Br7 equivalent widths for the 2 knots Kl and K2. Evolutionary tracks for starburst 
models with 3 different exponential decay timescales are also shown, reproduced from fig. 3 of Puxley, Doyon, & Ward (1997). The 
numbers on each track mark the starburst age in Myr. 
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